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STRESSES AT THE TIP OF A LONGITUDINAL CRACK IN A PLATE STRIP* 


By W. B. Fichter 
Langley Research Center 

SUMMARY 

The stress field near the tip of a central longitudinal crack in a plate strip is 
investigated. Several plane elastostatic problems involving opening and in-plane shearing 
of the crack, as well as one plate bending problem, are analyzed. Stress intensity fac- 
tors are obtained as a function of the ratio of strip width to crack length. In addition, for 
several problems some plots of the stresses on the longitudinal edges are presented to 
illustrate some additional effects of the proximity of the crack to the boundary. Finite 
strip width is shown to exert a strong influence on stress intensity factors in certain 
ranges of the ratio of strip width to crack length. Results for the fixed-edge plane prob- 
lems and the plate bending problem are only slightly influenced by changes in Poisson's 
ratio. The results should be useful in the design of fracture test specimens and in the 
analysis of fracture test data. 


INTRODUCTION 

In the aerospace industry there is concern over the deleterious effects of cracks 
in metal structures. Many in-service failures of aircraft components have been traced 
to fatigue cracks. In addition, the use of some high strength alloys has been compro- 
mised by crack sensitivity. Such considerations have motivated extensive development 
of the technology of predicting fracture of materials in the presence of cracks. An 
important aspect of this technology is the determination of stress fields near cracks 
and the study of the influence of loadings and boundary conditions on these stress fields. 

It is widely known that the linear theory of elasticity predicts a square root singu- 
larity in stress at the tip of a crack in a plate when the plate is subjected to loads which 
cause relative displacement of the opposing crack faces. (See, for example, refs. 1 to 3.) 
Although this singularity in stress is physically unrealistic, the strength of the stress 
singularity or some numerical modification of it, usually referred to as the stress 
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intensity factor, furnishes a means by which a few judiciously designed experiments 
might be utilized to predict the crack sensitivity of a larger class of configurations. As 
a result, requirements for costly and time-consuming testing programs might be 
reduced. 

Analytical information on stress fields around cracks in doubly infinite plates 
abounds in the technical literature. (See, for example, refs. 1 to 3, where numerous addi- 
tional references are cited.) On the other hand, relatively few analyses have been made 
of cracks in finite plates. More analytical information on the effects of finite in-plane 
plate dimensions and the range over which these effects are important would be useful in 
the design of practical structural components. 

Almost all analyses which have been performed on finite plates have been confined 
to cases in which the finite-plate dimension is parallel to the crack direction. The pres- 
ent report contains analyses of an infinitely long plate strip of width 2ah containing a 
central longitudinal crack of length 2a; that is, the finite-plate dimension is perpendic- 
ular to the crack direction. Several plane elastostatic problems and one bending problem 
are treated. Stress intensity factors are determined as a function of h, the ratio of 
strip width to crack length, to facilitate an assessment of the effects of finite strip width. 
The results are normalized with respect to those for a doubly infinite plate under the 
same loading conditions. For plane problems involving geometrical restraint of the 
longitudinal edges of the strip, some illustrative plots of stress distributions on the 
edges are presented. 

In all the problems analyzed herein, the external loading has been applied to the 
crack surfaces in order to make use of well-developed transform methods. However, 
with the exception of the problems involving concentrated forces, these problems can be 
easily extended by superposition to problems involving stress-free crack surfaces and 
uniform stresses or displacements on the longitudinal edges of the strip. The stress 
intensity factors are unaltered by such extensions. 

SYMBOLS 

A,B,C,D unknown coefficients in general solution of transformed biharmonic equation 
A m ,A n coefficients in series 

a half length of crack 

D plate bending stiffness 


2 


E 


Young's modulus of elasticity 


f(£) unknown function in dual integral equations 

G modulus of shear rigidity 

G(£) function in dual integral equations 

g(x) function of applied load in dual integral equations 

h nondimensional strip width 

K stress intensity factor 

k arbitrary real positive number in solution of dual integral equations 

L m,n coefficients in algebraic equations 

M nondimensional applied bending moment 

Mx>My>MxY nondimensional bending and twisting moments 

m,n integers 

P applied concentrated force 

Pm normalized series coefficients 

Qx>Qy nondimensional shear forces 

t plate thickness 

u,v nondimensional in-plane displacements 

w nondimensional lateral displacement 

x,y nondimensional rectangular Cartesian coordinates 

r gamma function 
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e € .v elastic strains 
x’y ,7 xy 

1 + JU. 


parameter in dual integral equations 




Poisson’s ratio 


transform variable 


nondimensional applied stress 


CT x ,Oy,a X y nondimensional stresses 


nondimensional applied shear stress 


4 > 

V 4 

Subscripts: 


nondimensional Airy stress function 


/ a 2 o 2\/ 8 2 a 2 \ 

biharmonic operator, — ^ ? — 5 + — s 

\9x 2 8y7\9x 2 8y 2 j 


refers to strip 


oo refers to doubly infinite plate 

Primes denote dimensional quantities. 

ANALYSIS 

The configuration is an infinitely long 
plate strip of width 2ah with a central 
longitudinal crack of length 2a. (See 
fig. 1.) The problems are divided into 
categories in accordance with the types of 
loading applied to the crack faces. In the 
first group of problems, the crack is being 
either opened or sheared by uniform loads. 


Bars denote transformed quantities. 


y' 



Figure 1.- Plate strip and coordinate system. 
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The second group involves the opening of the crack due to the application of concentrated 
loads to the crack faces. 


The plane problems and the bending problem are analyzed similarly. All the 
analyses are based on a method presented by Sneddon in reference 4 for determining the 
appropriate pair of dual integral equations to be solved for the Fourier transform of the 
Airy stress function. In the bending problem, the Fourier transform of the lateral dis- 
placement w is sought. After the dual integral equations have been obtained, a method 
given by Tranter in reference 5 is employed to solve the dual integral equations and 
eventually to obtain the stresses and displacements in the form of infinite series. Here, 
however, the primary goal is to obtain the stress intensity factors as functions of the 
ratio of strip width to crack length. 

\ 

For the plane problems the equations of plane strain are assumed to govern, it 
being understood that the plane strain solutions can be converted to plane stress solu- 
tions by the appropriate substitutions for Young's modulus E and Poisson's ratio j a. 
(See, for example, ref. 6.) 

Under the plane strain assumptions, the following system of equations governs the 
plane problems: 


where 


with the stress- strain relations 


vV( X ',y) = o 


8 2 <ft' 
9y' 


q2 j t 

<c-i4 - 


8x 


,2 



a 2 ft' 
ax' 9y_^ 


( 1 ) 


( 2 ) 



( 3 ) 
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and the strain-displacement relations 


t 



8u' 

ax’ 





av' 

9y' 





8y* 8u’ 

ax’ ay^ 


The quantities contained in equations (1) to (4) are made nondimensional by the 
following relationships: 

0* = EaV 

= Ea x 

Oy — E(7y 

<^xy = Eo xy 
6 x = e x 

€ y = € y > 

y xy = ^xy 
u’ = au 

v’ = av 
x’ = ax 
y’ = ay 

Equations (1) to (4) then assume the following form: 

v 4 0(x, y) = 0 


where 


,2 


with 


and 


€ X - °x “ M°y 
e y = CT y - m°x\ 


Vxy - “ a 


G ^ 


e y 


9u 

9x 

9 V 

3y 




_ 9v 9u 

y xy " a x + ay^ 


( 8 ) 


(9) 


In terms of the nondimensional coordinates, the strip is of width 2h and contains a 
crack of length 2. 


The fact that the plate is infinite in length strongly suggests the use of infinite 
integral transforms. With the complete Fourier transform of a function defined by (see, 
for example, ref. 4) 

f(£,y) = (* f(x,y) e -1 ^ x dx (10) 

^ - CO 

and inversely 

f(x,y) =^J oo f"a y ) e^ x d£ (11) 


equation (6) becomes 



<M4,y) = 0 


( 12 ) 


it being assumed that stresses and displacements vanish as |x| — The general solu- 
tion to equation (12) may be written as 


<?(£>y) = (A + By)e“^ y + (c + Dy)e^ y 


(13) 


where A, B, C, and D are functions of 4 to be determined by application of the 
boundary conditions. In terms of <f>, the transformed stresses and displacements are 
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and 











d£ 

d y> 





- i 


d 3 <j> 

dy 3 




( 14 ) 


(15) 


Because of similarities in the various analyses, two example problems are pre- 
sented in detail in the text, the remaining problems being discussed only briefly. 

Detailed analyses of the remaining problems are relegated to appendixes A and B. 

The first example problem concerns the opening of the crack by uniform pressure, 
the edges of the strip being taken to be free of stress. After this detailed treatment, 
brief discussion is given of other problems which can be analyzed in a similar manner 
such as uniform normal and shear stress distributions and uniform bending- moment dis- 
tribution on the crack faces. Detailed analyses of these other problems are given in 
appendix A. 

The second example problem concerns the opening of the crack by symmetrically 
applied concentrated forces, the edges of the strip being stress free. This problem is 
followed by brief discussions of other concentrated force problems, detailed analyses of 
which are given in appendix B. 

In the analyses, certain assumptions are made with regard to the permissibility of 
interchanging the order of summation and integration in some infinite series. These 
assumptions appear to be justified by the final results. 


Crack-Opening Mode — Uniform Pressure 

Free longitudinal edges .- In the first problem the edges of the strip are free of 
external load, and the crack is being opened by a uniform pressure Eo, where o is a 
nondimensional measure of the applied pressure. Because of symmetry, only the upper 
half of the strip (0 = y I h) needs to be considered. The boundary conditions in the 
dimensionless system are 
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a y (x,h) = 0 

°xy( x >h) = 0 

°xy( x >9) = 0 

a y (x,0) = -a 
v(x,0) = 0 

The transformed boundary conditions are 

Oy(£,h) = 0 
^xy(l>h) = 0 

®xy(£>0) = 0 




\ ( 16 ) 

(l*l<>) 

(M >l i 


■ a= iC^ feo)ellXd5 (m* 1 ) 

0 = j* ^ v(|,0) e 1 ^^ (|x|>l) 


In terms of the transformed stress function (p, the boundary conditions are 

feh) = 0 

^( { ,h)=0 

aft«.°) * » 

2iro = f ^a,0) e^ x dS 

*■'' — 00 

and, use of equation (19) having been made, 

0 = r&i,«)e^«| 

J -oo dy^ I 2 


(|x|<l) 

(M>») 


(17) 

(18) 

(19) 

(20a) 

(20b) 


Substitution of equation (13) into equations (17) to (19) yields 
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( 21 ) 


A = ^(e 2, ^ |h - 2| 2 h 2 + 2|||h - l) 

B = i^(e 2 ^ h + 2 m h-l) 

C =^(2| 2 h 2 + 2|?|h+ 1 - e' 2 ^l h ) 

where 

A = |4|(2|||h+1- e' 2 l^ h ) 


The direct stresses are certainly even functions of x and so, therefore, is the 
stress function <j>(x, y). It is then easily shown that 0(£,y) is an even function of 
This statement further implies that D is an even function of £, and allows equa- 
tions (20) to be modified to 

= f °° ^(cosh 2£h - 2£ 2 h 2 - l)cos |x d£ (| x | < l) (22a) 

2 J 0 A 

and 

0 = f -^(sinh 2£h + 2|h)cos £x d£ (|x|>l) (22b) 

J 0 A 


Equations (22) are the dual integral equations which must be solved for the unknown 
function D(|). 

In order to cast equations (22) in a form suitable for treatment by the method of 
reference 5, the following identity is used: 


cos £x = ^-l/2^ x ^ 

and the following substitutions are made: 

3/2 

f(|) = l-£-2(sinh 2£h + 2£h) 

> 

Q(t) = /cosh 2gh - 2£ 2 h 2 - l ) 

\ sinh 2£h + 2£h / j 


(23) 


(24) 


In equation (23), is the Bessel function of the first kind of order -1/2. With 

these substitutions, equations (22) assume the form 


g (x) = go-) m J_ l/ 2 (fr) (I x l< X ) < 25a) 
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and 


(|x | > l) (25b) 


where 



j _1/ 2 (£x) d S 



and the unknown function is now f(|). According to reference 5, equation (25b) is satis- 
fied automatically by 


f(l) = i Uk J A m J ,(f) (26) 

■SO 2m+k -2 

in which k is required only to be real and positive, and the A m values are to be 


- 1/2 

determined through the satisfaction of equation (25a). When g(x) is of the form Xx ' , 
further use of the method of reference 5 converts equation (25a) into the following infinite 
set of linear algebraic equations in A m : 


A 0 + l 

/< 

II 

a 

t 

•s 

a 

J 

m=0 


OO 
r — v 


A n + Z 

L m,n A m - 0 

m=0 



‘‘' Ms ) 

r H) 


/ (27) 


(n = 1, 2, 3, . . .) 


where 


L m ,n = ( 4n + 2k - 1) T [| 2 - 2k G(4) - l] J j(4) J ^ 

J 0 J 2m+k - i 2n+k-i « 

z z 


(28) 


For the present case of uniform pressure in the crack, X = (ir/2)^^a, and it is 

convenient to set k = 3/2 and A n = -^ P n . Equations (27) and (28) become 

z 


OO 

P 0 + I L m,O p m = 1 

m=0 

OO 

P n + ^ L m n P m = 0 (n = 1, 2, 3, . . .)| 


) (29) 


m=0 
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and 


L m,n = 2 < 2n + V J 0 - X ) J 2m+1 ( ^ J 2n+1 (?) ^ 
which for the present problem becomes 


^m,n 2(2n + 1). \ 


< 30 > 


0 


sinh 2ijh + 2|h 


For a specified value of h (>0), a square array of L mn values is computed from 
equation (30) by numerical techniques. The order of the array required for accurate 
determination of the P n terms (hence, D and eventually <#>(£, y)) is dependent on the 
assigned value of h. For large values of h (>1) the coefficients in equations (29) con- 
stitute a very nearly diagonal matrix. In fact, as h tends toward infinity, the matrix 
becomes diagonal, and the only nonzero root of equations (29) is P Q = 1/2. In general, 
then, the P n values decay more or less rapidly with increasing n, depending on whether 
h is large or small. Specification of a large value of h dictates the computation of 
only a few L m ,n terms, whereas a small value of h (<1) calls for a large array of 

L m , n terms. 

Once a sufficiently large array of L m>n terms has been computed, the L m>n 
terms are inserted into a truncated system of equations (29), which is then solved for Pn- 
In all the problems treated here, more simultaneous equations have been solved than 
appeared to be necessary for good convergence. The computations have been performed 
for numerous values of h in the range 0.1 = h ^ 6. 

On the line containing the crack (y = 0), 

cr y (x,0) = - £ 2 0(£,O) COS £jx d£ 

which by virtue of equation (13), equations (21), the first of equations (24), and equa- 
tion (26), becomes 


CT y (x,0) = -a 



cosh 2fh - 2£ 2 h 2 - 1 
sinh 2£h + 2i;h 


J 2m+1 ( ^ C0S * X d4 


or in a form better suited to present purposes 



2g 2 h 2 + 2^h + Li e — (g) cos £x d£ 
sinh 2£h + 2|h ) 2m+l 


( 31 ) 
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From reference 7, after slight simplification, 


f 00 cos[(2m + l)sin"^xl 

J 0 <*“»<■«- L 


A 

(x < 1) 


£ J 2m+1^ cos * x d * 


m+1 


> ( 32 ) 


( X 2 - l) 1 ^ + ( x 2 . i)l/2 


2m+l 


(x > 1) 


J 


The region of primary interest here is x > 1 which contains the area immediately ahead 
of the crack tip. From equations (31) and (32), 


CO 

<T y (x,0) = ct ^ P m 


il 11 


m 


(x2 - l) l / 2 [x + (x 2 - l) 1 / 2 ] 


It is apparent from physical considerations that the only singular point of oy(x,0) 
is at x = 1. It is now convenient to define a (dimensional) stress intensity factor by 


K’ = 2 1 / 2 lim f(x T - a) 1 / 2 o' (x\0)l 
x'-aL y J 


or in terms of dimensionless quantities 

K’ = Eo^a) 1 / 2 lim 
x-1 

where Ea is the uniform pressure applied to the crack faces, and the \/2 has been 
inserted strictly for convenience. Then the dimensionless stress intensity factor can 
be defined by 



K* 

Eoa 1 / 2 



lim f(x - 1) 1/2 J(x,0) 
x— 1L 


(34) 


It is shown in appendix C that the infinite integral in equation (33) is not singular at 
x = 1; therefore, it follows that the stress intensity factor is contained entirely in the 
first term on the right-hand side of equation (33) because only that term survives the 
limiting process indicated in equation (34). 
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Since the dimensionless stress 
intensity factor for an infinite plate 
under the same loading is given simply 
by 

Koo= 1 

then from equations ( 33 ) and ( 34 ) the 
ratio of dimensionless stress intensity 
factors for the strip and the infinite 
plate is 

y(-i)“p m os) 

m=0 

The ratio may be considered 

as a multiplying factor on the infinite- 
sheet stress intensity factor Koo to 
account for the effect of finite strip 
width. This ratio is plotted in 
figure 2. 



h 

Figure 2 .- Variation of Kp/Koo with h. Uniform pressure and free 
longitudinal edge. 


This problem can be changed, by 
superposition of a constant uniaxial 
stress field Oy(x,y) = ct, to a physically more realistic problem in which the crack sur- 
faces are free of applied load and the longitudinal edges y = ±h are subjected to a uni- 
form normal tensile stress. The stress intensity factor is unaffected by this alteration. 


Zero normal displacement of longitudinal edges .- This problem differs from the 
preceding one only in that the condition of zero normal stress on the longitudinal edges is 
replaced by the condition of zero normal displacement. With symmetry about the lon- 
gitudinal center line (y = 0) accounted for, the boundary conditions on the upper half of 
the strip are 


v(x,h) = 0 


CT xy (x,h) = 0 


°xy( x >0) = 0 


<< q 

3 

n 

Q 

(l*l <i) 

v(x,0) = 0 

(l*l >i) 
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This problem is solved by the same procedure employed in the first problem. The ana- 
lytical details are presented in appendix A. The ratio of stress intensity factors is again 
found to be 


K f 


K 


OO 


!«■*> 


m=0 


m. 


m 


Here, of course, the P m terms generally differ from those found in the first problem. 
The ratio Kp/Koo is plotted in figure 3 for 0.1 s h i 6. 

Another quantity of some interest is the normal stress distribution on the edge 
y = h. This stress arises from the zero normal displacement boundary condition. In 
appendix A, it is found to be 


c y (x,h) = -2<j l p m r Sinh c f sh %\ h h C _ 0Sh a J 2m+1 ({, cos ? x d ? 

m=0 * 


This stress is plotted in figure 4 for h = 2, 1, and 0.1 in the range 0lxi2. It is 
symmetric about the line x = 0. 


By superposition, this problem can be converted into one in which the crack faces 
are free of applied load and the longitudinal edges y = ±h are uniformly displaced away 
from the longitudinal center line; that is, in dimensionless terms, 

v(x,h) = -v(x,-h) = oh 


The dimensionless normal stress on 
the edges y = ±h required to pro- 
duce this uniform edge displacement 
is given by 

a y ,i = ct + Oy(x,h) 

This modified version, of course, 
yields the same stress intensity fac- 
tor and corresponds to constant- 
edge-displacement tensile tests if 
the testing machine grips do not 
introduce significant shear stresses 
at the edges. However, if the testing 
machine grips produce an effectively 
clamped- edge condition, the results 
of the following analysis should prove 
to be more pertinent. 


l 


.8 


.6 



.2 


0 


1 


± 

2 


_L 

3 


J 1 I 

4 5 6 



h 

Figure 3 .- Variation of Kp/Koo with h. Uniform pressure and zero 
normal displacement of longitudinal edges. 
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Clamped longitudinal edges . - 
For the cracked plate strip with 
clamped edges and uniform pressure 
in the crack, the boundary conditions, 
with symmetry accounted for, are 

u(x,h) = 0 
v(x,h) = 0 

<^jy(x,0) = 0 

o y (x,0) = -o (|x|<l) 
v(x,0) = 0 (|x|>l) 

By the procedures employed previ- 
ously, the ratio of stress intensity fac 
tors is again 



Figure 4.- Normal stress on longitudinal edge. Uniform pressure and 
zero normal displacement. 



h 

Figure 5.- Variation of Kf/Koo with h. Uniform pressure and 
clamped longitudinal edges. 


The P m terms differ from those of 
the preceding problems and, because 
of the clamped- edge condition, depend 
on the value assigned to Poisson's 
ratio ix. The ratio has been 

computed for fx = 0 and 0.5, two values 
which encompass those commonly used 
for most engineering materials, and 
for numerous values of h in the range 
0.1 s h s 6. This ratio is plotted in 
figure 5. 

It might be noted here that these 
results have been derived under plane 
strain assumptions. It is easily 
ascertained, however, that the results 
for plane stress are bracketed by the 
plane strain results presented in 
figure 5. 

The normal stress on a longi- 
tudinal edge is found to be 
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<V<X ’ h> = (T7 m>* 2 p m I" 2 C °3 h .^ (1 * ^ Sinh J 2mt i(9 i* d{ 

^ + m=0 0 — " ^ sinh 2£h - 2£h 

X *4“ 


This stress has been computed for ju = 0 and 0.5, and is plotted in figure 6 for 
h = 4, 1, and 0.1, in the interval 0 1x^2. 



Figure 6.- Normal stress on longitudinal edge. Uniform pressure and 
clamped longitudinal edges. 


Edge Sliding Mode - Uniform Shear 

The problems contained in this section are concerned with the sliding of the crack 
faces on each other in the plane of the plate due to the application of uniform shear stress 
Et to the crack faces. Again from symmetry considerations, it is possible to restrict 
attention to the upper half of the strip. 

Free longitudinal edges .- In this problem the longitudinal edges of the strip are 
free of external load. The boundary conditions on the upper half of the strip are 

cr y (x,h) = 0 

tfxyM) = 0 
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h 

Figure 7.- Variation of Kp^K^with h. Uniform shear stress and 
free longitudinal edges. 


Oy(X,0) = 0 

tfxyM) = T (M <:1 ) 

u(x,0) = 0 (|x|>l) 

In these shear problems, it is now the 
shear stress which is singular at the 
crack tip. Aside from this fact, how- 
ever, the analysis of this problem dif- 
fers little from those which precede it. 
The ratio of stress intensity factors is 
determined in appendix A to be 



This ratio of shear stress intensity fac- 
tors has been calculated for 0.1 S h ^ 6 
and is presented in figure 7. 


Zero normal displacement of longitudinal edges. - This problem differs from the 
preceding one in that the longitudinal edges of the strip are restrained against normal 
displacement. The boundary conditions on the upper half of the strip are 

v(x,h) = 0 

°xy( x > h ) = 0 
Oy(X,0) = 0 


°x y( x >°) = T (l x l <:l ) 

u(x,0) = 0 (| X I >:1 ) 


Details of this analysis also are contained in appendix A. The ratio of stress intensity 
factors has the same form as in the preceding problems. It has been computed for 
0.1 and is presented in figure 8. 

Restriction of normal displacement of the longitudinal edges gives rise to normal 
stress there. This stress is found to be 


Oy(x,h) = r 



£h sinh gh 
cosh 2£h + 1 


J 2m+1 ( ^ sin 
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This normal stress is plotted for 
h = 2, 1, and 0.1 in figure 9, in the 
range 0 ^ x ^ 2. It is an antisym- 
metric function of x. 

The stress intensity factor for 
this problem is identical to the stress 
intensity factor for the problem in 
which the crack faces are free of 
stress, and the longitudinal edges are 
subjected to uniform shear stress 
while being restrained against normal 
displacement. 

Zero tangential displacement of 
longitudinal edges .- In this problem the 
longitudinal edges of the strip are 
restrained against tangential displace- 
ment and the crack faces are subjected 
to uniform shear stress Et. The 
boundary conditions on the upper half 
of the strip are 

<r y (x,h) = 0 

u(x,h) = 0 

oy(x,0) = 0 

°xy( x >°) = T (|x|<l) 

u (x,0) = 0 (| x | > l) 

For this problem the ratio of stress 
intensity factors again has the same 
form. It has been computed for numer- 
ous values of h in the range 
0.1 ^ h = 6. The results are presented 
in figure 10. 

Shear stress is present at the lon- 
gitudinal edges of the strip, because of 
restriction of tangential displacement. 

In appendix A, this stress is found to be 



h 

Figure 8.- Variation of Kp/K^ with h. Uniform shear stress and 
zero normal displacement. 


y 


i T pbrH 



Figure 9.- Normal stress on longitudinal edge. Uniform shear 
stress and zero normal displacement. 
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V*>h) - * l Pm ***£&'!? * W» cos * d 5 

m=0 


u = o y =0 


It has been computed for h*2, 1, 
and 0.1 and is plotted in figure 11 for 
0 ^ x g 2. It is symmetric in x. 

The stress intensity factor for 
this problem is the same as for the 
problem in which the crack faces are 
stress-free and the longitudinal edges 
are uniformly displaced tangentially by 
an amount u(x,±h) = ± ■ and are 

(jT 

free of normal stress. 

Bending Problem 

The lateral bending of a thin 
strip containing a central longitudinal 
crack gives rise to infinite moments 
(hence, infinite bending stresses) at 
the tips of the crack, according to the 
classical small- deflection theory of 

plates. In dimensionless form, the equations governing bending of a thin plate in the 
absence of lateral load are (see ref. 8) 

V 4 w(x,y) = 0 



Figure 10.- Variation of Kp/K^ with h. Uniform shear stress and 
zero lateral displacement. 


Moments: 


M X = 


My = 



+ M 


+ M 


Shear forces: 


M XY = " M YX = ^ 


Qx = -lfe v2w 

Qy = • Jb v2 * 


20 






X 

Figure 11.- Shear stress on longitudinal edge. Uniform shear stress 
and zero lateral displacement. 


In the present problem, the longitudinal edges of the strip are simply supported and sub- 
jected to a uniform bending moment, designated as -M in the dimensionless system. 
For purposes of determining the stress intensity factor, a uniform bending moment 
M Y (x,y) = M is superimposed on the original system; this process results in uniform 
bending moments on the crack surfaces and longitudinal edges which are free of bending 
moment in the y-direction. Attention again is restricted to the upper half of the strip 
because of symmetry considerations. 


The boundary conditions are 


w(x,h) = 0 


My(x,h) = 0 


_ 9 _ 

9y 


%x,0) ♦ (2 - v) 

|_9y 2 




= o 


My(x,0) = M (|x| < l) 

^ x ,o) = o (|x|>l) 
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Because the bending problem is also governed by the biharmonic equation, the method of 
analysis employed in the preceding problems is applicable. The quantity of primary 
interest is now the bending moment My(x,0). In appendix A, it is found to be 


00 

M y (x, 0) = -M ^ P m 


Ml 


m 


„feo k2- 1 )'/2[ x + (x 2. l) l/2]— 

r °o 1 - 2|h(i-=-^+ e" 2 £ h 

+ \ U j 0 .(^cos^xdl 

Jn cosh 2£h + 1 2m+1 


jq cosn 2£h 

The ratio of bending stress intensity factors for the strip and the infinite plate then is 


K f 

K~ 


= £ (-D m p 


m 


m=0 


Note the dependence of My(x,0) on Poisson's ratio /i. The ratio K^/K^ has been 
computed for ja = 0 and 0.5 for numerous values of h in the range 0.1 = h = 6. The 
results are presented in figure 12. 



Figure 12 .- Variation of Kp/Kc© with h for bending problem. 


Crack Opening Mode - 
Concentrated Forces 

In the three problems which follow, 
the crack faces are subjected to equal and 
opposite compressive concentrated forces 
at their centers. Because the analysis of 
this particular class of problems differs 
somewhat from preceding analyses, the 
first of these problems is analyzed in 
some detail, and the analysis of the 
remaining problems is relegated to 
appendix B. 

Free longitudinal edges .- In this 
problem the edges of the strip are free of 
external load and the crack faces are sub- 
jected to concentrated forces of (nondi- 
mensional) intensity P/Ea. Attention is 
restricted here to the upper half of the 
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strip by virtue of symmetry properties. The boundary conditions on the upper half of 
the strip are 

CTy(X,h) = 0 

°xy( x >h) = 0 
cr^foO) = 0 


CT y (x ’ 0) = "Is 6( x) (l x l<i) 

v(x,0) = 0 (| x | > l) 


where 6(x) is the Dirac delta function. The transformed boundary conditions are 

?(5,h) =0 


^<{,h) = 0 

^ 5 ,°) = 0 

= j f { 2 3>({,0)cos {x d£ 


> (36) 


(M<») 



(l*M 

J 


Substitution of the general solution for <p (eq. (13)) into the first three of equations (36) 
again yields equations (21). The following substitutions are made 

cos |x = / J_ 1 y 2 (|x) 

3/2 

f(£) = - - (sinh 2£h + 2|h) 

G(t) = t /cosh 2£h - 2£ 2 h 2 - l \ 

\ sinh 2|h + 2£h / 


The dual integral equations then assume the form 

g(x)=J G(«)f(«) J_ 1/2 (£x) d| 


(| x | < l) (37a) 
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and 


0= Ct(& J.i/2^ x ) d ^ 


(|x| > l) (37b) 


r- 1/2 


where g(x) now has the form Ax’ ' 6(x), Once more it is assumed that 


f(!) - S 1//2 \ A m J 2m+1^ 
m=0 

which automatically satisfies equation (37b). Further application of the method of refer- 
ence 5 yields a different set of linear algebraic equations 

OO 

Pn + I Wm = 1 <" ' 2 > < 38 > 

m=0 


where 


and 


F P — A 
Ea Pn An 


L m,n = 2 < 2n + ^ J n “ *) J 2m+1^ J 2n+1^ 


V * 

For x > 1, the stress Oy(x,0) is found to be 


°y (x ’ 0)_ |fe I Pm 


m=0 


f 

(-I) 1 " __ 

l(x 2 - l) 1/2 

x + (x2-l) 1/2 


2m+l 


+ f 00 2 ^ 2 h 2 + ± 1 - e- . - . v . J (|) cos |x d|) 

Jo sinh 2£h + 2£h 2m+1 < 


»-2|h 


(39) 


For the problem of concentrated forces it is again convenient to express the effect 
of strip width in terms of the ratio of stress intensity factors. For the infinite plate 
with the same loads, the stress o ^foO) is given by 


'y>°° 

oo 


O y>OO (X,0) = 2 ^ a £ 




m 


(x 2 - l) 1/2 |x + (x 2 - lH 


,a/2l 


2m+l 


(40) 
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Then the ratio of stress intensity factors 
may be written as 


K f o- (x,0) 

— - = lim — - — - — - 
Kco x-1 a oo(x,0) 


(41) 


After substitution of equations (39) and (40) 
into equation (41) and appropriate passage 
to the limit, the ratio of stress intensity 
factors is found to be 


g=l + 2 2 (-l) m (P m - 1) (42) 

m=0 

This ratio has been computed for 
0.1 ^ h i 6, and the results are plotted 
in figure 13. 

Zero normal displacement of 
longitudinal edges .- This problem differs 
from the preceding one in that the longi- 
tudinal edges are restrained against nor- 
mal displacement. The boundary condi- 
tions on the upper half of the strip are 



h 


Figure 13.- Variation of Kp/Koo with h. Concentrated force 
and free longitudinal edges. 


v(x,h) = 0 
o-xyfoh) = 0 


CT xy( x >0) — ® 


a y (x,0) = - ^ 6(x) (|x| < l) 

v(x,0) = 0 (|x | > l) 

The ratio of stress intensity factors for this problem is found in appendix B to have the 
form of equation (42). The results are presented in figure 14. 

The normal stress on the edge y = h also has been obtained. It is 


"AW - - S l P» f Sinh c f sh ^rf ?h W® cos IX dj 


m=0 


This stress is plotted in figure 15 for h = 4, 2, and 1. 
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Clamped longitudinal edges .- In 
this problem the longitudinal edges 
are completely restrained against 
normal and tangential displacements. 

The boundary conditions on the upper 
half of the strip are 

u(x,h) = 0 
v(x,h) = 0 

°xy( x >0) = 0 

o y (x,0) = - ^ 6(x) (| x | < l) 

v(x,0) = 0 (|x|>l) 

The analytical details of this problem 

are contained in appendix B. The Figure 16.- Variation of Kp/Koo With h. Concentrated force and 

ratio of stress intensity factors again clamped longitudinal edges. 

has the form of equation (42). It has 

been computed for numerous values of h in the range 0.1 i h ^ 6, and the results are 
presented in figure 16. For this problem, the ratio of stress intensity factors depends 
on Poisson’s ratio. It has been computed for \i = 0 and ji = 0.5. 

There are stresses along the longitudinal edges because of the clamped-edge con- 
dition. The normal stress on an edge is found in appendix B to be 



cr (x,h) = ^P 5 Y P m f 2 cosh gLtlLt sin! 

tfEa(l + n) 2 £ ^0 sinh 2 _ 2 

1 + JLt 


J 2m+1^ cos 


This stress is plotted for h = 4, 2, and 1 in figure 17 for \± - 0 and p = 0.5. 


RESULTS AND DISCUSSION 


It is convenient to consider the ratio K-p/Koo as a multiplying factor on the 
infinite- sheet stress intensity factor to account for the finite width of a plate in an actual 
structure or to account for the finite distance between testing machine grips in laboratory 
tests of cracked sheet. With regard to the latter application, the results presented 
should prove useful in the correlation of existing experimental data. These results also 
should be helpful in the design of test specimens, because they furnish the specimen 
designer with guidelines for minimizing the perturbing influence on the stress intensity 
factor of the boundary conditions imposed by the testing machine grips. 
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The Kp/K^ curves for the six 
crack-opening mode problems are pre- 
sented in figures 2, 3, 5, 13, 14, and 16. 

For the problems in which the longitudinal 
edges of the strip are stress-free (see 
figs. 2 and 13), the ratio of stress inten- 
sity factors is seen to be always greater 
than 1, and to approach 1 rather rapidly 
for moderately large values of h, the 
ratio of strip width to crack length. On 
the other hand, the ratio is always less 
than 1 for the problems involving restraint 
of the normal, or normal and tangential, 
displacements of the edges of the strip. 

(See figs. 3, 5, 14, and 16.) Again, the 
ratio is seen to approach 1, although 
somewhat less rapidly than for free edges, 
as h becomes large. 

Figures 5 and 16 are for problems involving completely fixed longitudinal edges, a 
set of boundary conditions which causes the stress intensity factor to depend on Poisson's 
ratio. For these problems, the stress intensity factors have been computed for tx = 0 
and fj. = 0.5, two values of Poisson’s ratio which encompass the values commonly 
assigned to most engineering materials. The results suggest that Poisson's ratio exerts 
little influence on stress intensity factors for most practical crack-opening problems. 

In the crack- opening mode problems involving restraint of displacements at the 
longitudinal edges of the strip, normal stresses arise there. These normal stresses have 
been computed for various values of h and are presented in figures 4, 6, 15, and 17. 
Figures 4 and 6 are for problems involving uniform pressure in the crack, and figures 15 
and 17 are for the concentrated force problems. For larger values of h the normal 
stresses are rather small and well-behaved, but these stresses and their gradients grow 
as h becomes small because the stresses due to the applied forces and the crack-tip 
singularity have less area in which to diffuse and decay. 

Figures 6 and 17 are for the problems in which the longitudinal edges of the strip 
are fixed so that the stresses depend on Poisson's ratio. Much like the stress intensity 
factors for these problems, the normal stresses are not strongly influenced by changes 
in Poisson’s ratio. 

A curious result coming out of the two concentrated force problems which involved 
restraint of normal displacements at the strip edges can be seen in figures 14 and 16. 


0y(x,h) 

TiPT 

~ 7iEa 



Figure 17.- Normal stress on longitudinal edge. Concentrated force 
and clamped longitudinal edges. 
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As h becomes small, the stress intensity factor goes to zero at a seemingly arbitrary 
value of h, specifically, in the neighborhood of h = 0.8, the precise crossing varying 
with Poisson's ratio in the case of fixed longitudinal edges. This behavior is due to 
closure of the crack tips. The Kp/K^ curves have been terminated at the crossings 
because for smaller values of h the problems cease to be only crack problems, and 
instead become very complicated contact problems which are not expected to yield addi- 
tional significant information on the stress intensity factors. 

The Kp/Koo curves for the edge- sliding mode (shear) problems are presented in 
figures 7, 8, and 10. For the plate with free longitudinal edges (fig. 7), the stress inten- 
sity factor increases without limit as h approaches zero, and approaches the infinite- 
sheet stress intensity factor as h becomes large. This result is similar to those 
obtained for the crack-opening problems with the same free- edge conditions. 

For the second shear problem, in which normal displacement of the longitudinal 
edges is prevented, the stress intensity factor for the strip (fig. 8) differs little from its 
infinite- sheet counterpart except for very small values of h, where the problem has lost 
most of its physical importance. The restriction of normal displacement of the longi- 
tudinal edges gives rise to normal stresses there. This normal stress is presented in 
figure 9 for several values of h. Not unexpectedly, the normal stress is more severe 
for small values of h. 

An unusual result from the third shear problem, in which the tangential displace- 
ment of the longitudinal edges is prevented, is apparent in figure 10. The ratio Kp/Koo 
approaches 1 at both ends of the h-axis, and is greater than 1 at all interior points. This 
result means that there are two values of h corresponding to every admissible value 
of Kp/Koo. 

In this shear problem, shear stresses develop at the longitudinal edges. This 
shear stress is presented in figure 11 for several values of h. As in the preceding 
problems, the strip width exerts a strong influence on the stress at the boundary. 

For the one plate bending problem treated here, the Kp/K^ curves are presented 
in figure 12, where it is seen that the ratio of moment (or outer-fiber bending stress) 
intensity factors depends on Poisson's ratio. Again, the effect of Poisson’s ratio on 
Kp/K OO is small. The ratio of stress intensity factors approaches 1 for large h and 
increases without limit as h approached zero. Finite strip width is seen to have little 
influence for ratios of strip width to crack length greater than about 2. 

For the problems treated herein, finite strip width is seen to be of little significance 
for ratios of strip width to crack length greater than 6. In fact, for problems involving 
free longitudinal edges, h greater than about 3 is large enough to allow treatment of the 
strip as an infinite plate, at least for the purposes of crack-tip stress field investigations. 
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CONCLUDING REMARKS 


Stress intensity factors have been obtained for several problems involving crack- 
opening or edge sliding modes and, in one case, opening of the crack by bending, for a 
plate strip containing a central longitudinal crack. The stress intensity factors, normal- 
ized with respect to their counterparts for an infinite plate under the same loading, have 
been presented as functions of the ratio of strip width to crack length for the entire range 
of practical interest. For problems in which the stresses depend on Poisson's ratio, the 
results have been obtained for at least two values of Poisson's ratio which encompass the 
range of primary engineering interest. 

In addition, where restrictions of displacements of the longitudinal edges of the 
strip have given rise to boundary stresses, some sample stress distributions have been 
presented to illustrate their dependence on the ratio of strip width to crack length. 

For the problems treated herein finite strip width is seen to be of little significance 
for ratios of strip width to crack length greater than 6. In fact, for problems involving 
free longitudinal edges, h greater than about 3 is large enough to allow treatment of the 
strip as an infinite plate, at least for the purposes of crack-tip stress field investigations. 

The results presented herein should be useful in the conversion of data from tests 
of cracked sheet- metal specimens to practical information suitable for use by design 
engineers. The results should also be of assistance in the design of cracked- sheet test 
specimens to avoid edge effects. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., February 1, 1967, 

124-08-06-07-23. 
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APPENDIX A 


PROBLEMS INVOLVING UNIFORM LOADS ON CRACK FACES 

With the exception of the first crack-opening problem, which has been given detailed 
treatment in the body of this report, this appendix contains the bulk of the analytical 
details of the problems involving uniform loading of the crack faces. The three types of 
uniform loading considered are: (1) normal stress, (2) shear stress, and (3) bending 
moment. 


Crack Opening Mode — Uniform Pressure 

Zero normal displacement of longitudinal edges .- In this problem, the edges of the 
strip are restrained against normal displacement and the crack faces are being separated 
by uniform normal stress. Because of symmetry considerations, only the upper half of 
the strip needs to be considered. The boundary conditions there are 

v(x,h) = 0 

Qgy(x,h) = 0 
°xy( x ’0) = 0 


Oy(x,0) = -a 
v(x,0) = 0 

In terms of $(£,y) the transformed boundary conditions are 

^|(l, h) = 0 

dy 3 

* 0 

$ 1 , 0) -0 


(1*1 < l) 
( 1*1 > 1 ) 




>(A1) 


o = ^ £ 2 4>(£,0) cos £x d£ (M<l) 

0 = f (£,0) cos £x (| x | > l) 

J 0 dy 3 ¥ v j 
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Substitution of the expression for 4>(%,y) (eq. (13)) into the first three of equations (Al) 
yields ^ 

A = - ^(e 2 l*l h + 2|$|h - l) 


where 


B = -l£!2(e 2 l«l h -l) ) 

C = - 2(2|{|h + 1 - e' 2 l ? l h ) 
A = |t|(l - e- 2 l«l h ) 


Then by use of equations (A2), the identity (23), and the following substitutions: 

3/2 

f(£) = - £ -—-(cosh 2£h - 1) 

A 

G (t) = J sinh 2|h+ 2&i 
cosh 2|h - 1 j 

the last two of equations (Al) assume the form 

g(x) = J f(£) G(£) J _ 1 / 2 (?x) d| (|x| < l) 

°=J o f(l) J. 1 / 2 (? x ) (| x |>l) 

and again 

g(x) = ^ x" 1 ^ 2 

Therefore, the method of reference 5 may be applied. It is assumed that 

00 

Z PmJ 2m+l^ 

m=0 


(A2) 


} (A3) 


and the following infinite system of linear algebraic equations results: 
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Pn + ) oPm = * 


m=0 

OO 


*n + ^ ® (n = 1, 2, 3, . . .) 


> (A4) 


m=0 


L m,n " 2 ( 2n + D J Q - *) J 2m+1 ( ^ J 2n+1 ( ^ ^ 

After determination of the P n terms the stress Oy(x,0) can be written 


oo 

a y (x,°) = a ^ P m 


ill) 


m 


m=0 


xl/2f 


.4/21 


( X 2-ir / \x + (x 2 -l) 1 ^J 


2m+l 


-I 




(A5) 


As before, the infinite integral in equation (A5) can be proved convergent at x = 1. Con- 
sequently, the ratio of stress intensity factors is again found to be 

K f 


S' I ( - 1)mp - 


m=0 


Because of the restriction against normal displacements of the longitudinal edges, 
there is normal stress on these edges. In terms of the transformed stress function <$>, 
this normal stress is 


<*y(x,h) = - J |2 0(£ >h ) cos | X 

which with the proper substitutions becomes 

V x ’ h) - - 2 » 1 p m [ J 2m + 1® c ° s < A6 > 

m=0 

This stress distribution is symmetric in x. With a sufficient number of the P m values 
already computed for this problem, the computation of Oy(x,h) is merely a matter of 
numerically performing the integration indicated in equation (A6). 
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By superposition, this problem can be converted into a problem in which the crack 
faces are stress-free and the longitudinal edges of the strip are uniformly displaced 
away from the longitudinal center line; that is, in dimensionless terms, 

v(x,±h) = ±oh 

The (dimensionless) normal stress on the edges y = ±h which is required to produce 
this uniform edge displacement is given by 

°y, 1 ( x , ± h) = u+ °y( x > h ) 

This modified version, of course, yields the same stress intensity factor and corresponds 
to constant-edge-displacement tensile testing machine experiments if the testing machine 
grips do not introduce significant shear stresses at the edges. On the other hand, if the 
testing machine grips do produce a truly clamped-edge condition, the results of the fol- 
lowing analysis should prove to be more pertinent. 

Clamped longitudinal edges .- For the centrally cracked plate with clamped edges 
and uniform pressure in the crack, the boundary conditions for the upper half of the strip 
are 

u(x,h) = 0 
v(x,h) = 0 

°xy( x >°) = 0 

o y (x,0) = -a (|x|<l) 

v(x,0) = 0 (|x|>l) 


In terms of <p, the transformed boundary conditions are 

+ M£ 2 £(4,h) = 0 

dy^ 

^U,h) - (2 + m)£ 2 ®$,h) = 0 

dy” ^ 

44({,o) = o 

{ 2 feO)cos{xd5 

17 J o 

0 = f cos £x d^ 

J 0 dy J 



34 



APPENDIX A 


Substitution of equation (13) into the first three of equations (A7) gives 


where 





J 



(A8) 


and 


1 + JU 


By use of equations (A8), the identity (23), and the following substitutions: 

fc3/2 T 


m = - sinh 2 * h - 


/?-=-£ cosh 2£h + 2£ 2 h 2 + -i*- + A 
G(|) = | LLR L±J£ 


3 - jj. 

1 + li 


sinh 2£h - 2£h 


/ 


the last two of equations (A7) again take on the form of equations (A3). Application of the 
method of reference 5 yields the now-familiar system of equations (A4), which are to be 
solved for the P n terms. 


After some substitutions and manipulations, the stress Oy(x,0) for x > 1 is 
found to be 


OO 

°y( x ,0) = ® ^ P nil 


ill! 


m 


m=0 1 (x 2 - l) 1/2 [x + (x 2 - l) 1/2 ~ 


2m+l 



1 + + 2|h + 2£ 2 h 2 + e' 2 ^ h 

1 + M 1 + /1 

sinh 2£h - 2£h 

1 + ii 


J 2m+1^ cos 
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From this equation, it is found that the ratio of stress intensity factors has the form 



m 


An additional complication in this problem is that the stress intensity factor now 
depends on Poisson's ratio ix. In order to obtain an envelope of Kp/K^ curves, the 
L m , n terms (hence, the P m terms and Kp/K^) have been computed for /x = 0 and 
fx = 0.5, two values which encompass those commonly used for engineering materials. 

Because of the dependence on Poisson's ratio, the results for plane stress condi- 
tions will differ slightly from the present results, which are for plane strain. It is 
easily shown, however, that the Kp/K^ curves presented here bracket the curves for 
plane stress. In view of the generally small differences between the fx = 0 and fx = 0.5 
curves, it is not felt that the additional calculation of corresponding curves for plane 
stress is necessary. 

Again there is normal stress on the edges y = ±h. It is found to be 


°y( x > h ) = - 


4a 

(1 + M ) 2 



2 cosh $h + (1 + ju)$h sinh £h 
sinh 2£h - 2|h 

1 + ix 


J 2m+1^ cos d £ 


The required operations are performed numerically. 


Edge Sliding Mode - Uniform Shear 


The problems contained in this section are concerned with the sliding of the crack 
faces on each other in the plane of the strip due to the application of opposing uniform 
shear stresses Er to the crack faces. Again, because of symmetry considerations, 
the investigation can be restricted to the upper half of the strip. 

Free longitudinal edges .- In this problem the edges of the strip are stress free. 
The boundary conditions on the upper half of the strip are 

Oy(x,h) = 0 
o*y(x,h) = 0 
Oy(X,0) = 0 


°xy( x >0) = t ( I x | < l) 

u(x,0) = 0 (|x|>l) 
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The transformed boundary conditions are 

5(4, W = o 
f (5,h) = o 

= o 

T = ‘il 0 4 ^(4,0) cos 4x d{ (|x|<l) 

°= y o ^|(4,°) cos 4 X f (Ix|>l) 


Substitution of equation (13) into the first three of equations (A9) gives 

A = -C = - 2ll]h?D 

A \ 

B = - g ( e 2 N l h _ 2|^|h - l) 

J 

where 


A = 1 - 2|£|h - e' 2 l^l h 

With the use of equation (23), equations (A10), and the substitutions 

1/2 

f(£) = - H ... D .( S inh 2£h - 2£h) 

A 

<*> ■ *) 


(A9) 


(A10) 


the last two of equations (A9) assume the form of equations (A3), where 

*<«> ■ < f ) 1/2 ^ 1/2 

The method of reference 5 is now applied, and again equations of the form (A4) are 
obtained. 

However, for these shear problems, it is now the shear stress a xy (x,0) which is 
singular at the crack tip. For x > 1, it is found to be 
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CO 

a x y (x ’0) = -T Z P m 


i=n 


m 


m=0 l(x 2 - + (x 2 - l) 1/2 ] 


2m+l 


C" 1 - 2|h + 2 A 2 - e' 21 * , , .] 

J 0 slnh 2|h - 20. 


With the following definition for the shear- stress intensity factor 


K = 2 1 / 2 lim 
x—1 


(x- 1) 


1/2 


^(x,0) 


the ratio of stress intensity factors is 





m=0 


m. 


m 


The problem is equivalent, insofar as stress intensity factors are concerned, to one in 
which the crack faces are stress free and the longitudinal edges are subjected to uniform 
shear stress of magnitude Et. 

Zero normal displacement of longitudinal edges .- This problem differs from the 
preceding one in that the longitudinal edges are restrained against vertical displacement. 
The boundary conditions are 


o 

II 

xT 

s 


$ 

V 

n 

o 


Oy(X,0) = 0 


°xy( x > 0 ) = T 

(M <>) 

u(x, 0 ) = 0 

A 

X 
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The transformed boundary conditions are 

= 0 

dy J 

$<«.!>> - » 

0(€,o) = o 

«^te0)cos|xdj 

° = ^ a,0) cos {x f 




} (AH) 


(1*1 <i) 
(l *i > i) 


Substitution of equation (13) into the first three of equations (All) yields 

A = -c = M 

B = - E(e 2 l«l h + l)J 

where 

d = l + e- 2 l«l h 

With the use of the identity (23), equations (A12), and the following substitutions: 

f(£) = 2i | 1/2D (cosh 2£h + 1) 

A 

G (t\ = J sinh 2jh + 2gh N 
v?; A cosh 2£h + 1 


(A12) 


the last two of equations (All) take on the form of equations (A3). The method of refer- 
ence 5 is applied and equations (A4) again are obtained. For x > 1, the shear stress 
o xy (x,0) is obtained as 
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OO 

7 xy ( X ,0) =-T £ 


illl 


m 


m=0 j( X 2- l) l/2 [x + (x2. 1) 1/2 ] 


2m+l 


+ 



1 - 2gh + e~ 2 ^ h j (£) 
cosh 2£h + 1 2m+l 


cos £x d£ 


From this equation, the ratio of stress intensity factors is found to have the same form 
as in the preceding problems. 


The normal stress on the longitudinal edges, caused by the v(x,±h) * 0 boundary 
condition, is found to be 


°y( x ,h) = T f P m P J 2m+ 1^) Sin & 

y ^ cosh 2£h + 1 iim+i 

m=0 


In this case, the normal stress is an antisymmetric function of x. 

With respect to the stress intensity factor, this problem is equivalent to one in 
which the crack faces are stress free, and the edges of the strip are subjected to uniform 
shear stress while being restrained against normal displacement. 

Zero tangential displacement of longitudinal edges .- In this problem the edges of 
the strip are free of normal stress and restrained against tangential displacements. 

The boundary conditions are 

°y( x ,h) = 0 

u(x,h) = 0 

Oy(x,0) = 0 


°xy( x >°) = T (l x | < l) 

u(x,0) = 0 (|x| > l) 


In terms of 4>, the transformed boundary conditions are 

_ ~\ > 

<P(Z, h) = 0 

^f(£,h) =0) 

dy z j 

fe 0)=0 J ^ 

(Equations continued on next page) 
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T= ‘il 0 S ^ { ’°> cos d « 

> 

(l x l< l) 

0=f”^l({,0)cos{x^ 
J 0 dy 2 £ 

(l x l > l) 
J 


(A13) 


Equation (13) substituted into the first three of equations (A13) yields 


with 



(A14) 


A = l-e- 2 ^ h 

After substitution of equations (A14) into the last two of equations (A13), use of the iden- 
tity (23), and use of the following substitutions: 

f (£) = - ~- A D (cosh 2£h - 1) 


G(m = t fsinh 2gh - 2|h\ 

? \ cosh 2$h - 1 ) 

the system of equations (A3) again is obtained. The method of reference 5 is employed 
to solve these equations for the P m values, and the stress a X y(x,0) is 


OO 

°xy( x ,0) = -t £ P m 


illl 


m 


m=0 |(x 2 - 1) 1/2 ^ + (x 2 - i) 1//2 ] 


2m+l 



1 - 2gh - e~ 2 ^ h 
cosh 2£h - 1 


J 2m+1 ( ^ cos 


The ratio of stress intensity factors is again found to be 



m. 


m 
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The shear stress at the edge y = h, caused by restriction of the tangential dis- 
placement there, is given by 

V (x ’ h) ■ 2t 2 p ” C Sinh cth !{h TT jh «» 

m=0 

This shear stress is symmetric in x. The stress intensity factor for this problem is 

the same as for the problem in which the crack faces are stress free, and the edges of 

Exh 

the strip are uniformly displaced tangentially by an amount u(x,±h) = ± and are 
free of normal stress. 


Bending Problem 


This problem is concerned with the lateral bending of a thin plate strip containing 
a central longitudinal crack. In the classical small- deflection theory of plates (see 
ref. 8), the governing equations for a plate free of lateral load are 

DV^w'fr'jy’) = 0 


where 



and w' is the lateral deflection of the plate middle surface, with moments given by 




+ M 






m xy “ " M YX " ^ “ 


n) 


9 2 W' 
8x’ 9y' 


and shear forces given by 

«k - - 5 W V ' 2v ' 

< 5y = - 5 #' v ' 2w ' 

The foregoing equations can be nondimensionalized conveniently in the following manner. 
Let 
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w’ = aw 
x' = ax 
y’ =ay 

M k = J M X 


M Y “ a M Y 


m xy “ ~ M YX " a m xy = ' r m yx 


Q k=^x 


Q Y “ 2 q y 


Then the governing equations become 


where 


V w(x,y) = 0 


V 4 = a 4 V' 4 


with moments 


M X = 


My = 


^ 2 w + n 9 2 w\ 


\.dx 


9y 


L2 a 2 \ 

9 W 9 W) 

fJ. 


\9y 


9x" 




a^w 


M XY " " M YX ^9x 9yJ 


and shear forces 


Q x = -|; v2 W 


q y .-±v\ 


(A 15) 


(A 16) 


(A 17) 


In the problem to be analyzed here, the longitudinal edges of the strip are simply sup- 
ported, and a uniform bending moment DM/a is being applied to the crack surfaces 
(M is a dimensionless measure of the magnitude of the applied bending moment) . 
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In terms of dimensionless quantities, the boundary conditions are 


w(x,±h) = 0 
My(x,±h) = 0 


a 

3^W/ 

W 

lv‘ 


)^(x,0) 

dx* 


= 0 


M y (x,0) = M (l X | < l) 

^(x,°)=0 (|x|>l) 


The first two boundary conditions describe the simple supports. The third derives for 
x < 1 from the fact that the crack faces are free of shear force, and for x > 1 from 
considerations of symmetry. The fourth is due to the applied uniform bending moment 
(x < 1) and to symmetry (x > 1). 

After Fourier transformation, the governing equations (A15) to (A17) take the form 

(a? ' |2 ) 5<{,y) = 0 <A18) 


with moments 




M 


XY 


= -M YX =-i{(i- M >;| 


(A19) 


and shear forces 



(A20) 


In terms of w(£,y), the transformed boundary conditions, with symmetry about the lon- 
gitudinal center line accounted for, are 

w(£,h) = 0 (A21) 
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^■^ 2 K h) = 0 


(A22) 


d_ 

dy 


~~n - (2 - m)4 2 |w(^, 0) = 0 
dy^ 


(A23) 


-M7T 


-r 

Jn 


2 

4w(J,0) - M| 2 w(^,0)| 

dy^ 1 


cos |x d£ (| x | < l) (A24a) 


and 


0 = l~af (? ’ 0)cos5xd{ 


d y 

The general solution of equation (A18) can be written as 

w(£,y) = (A + By)e"1 ^ y + (C + Dy)e^ ^ y 
Substitution of equation (A25) into equations (A21) to (A23) gives 

A = §(l + 2|{|h + e 2 l^ h ) 

B = -li^(l + e 2 l«l h ) > 
C=-g(x-2| 5 |hK + e- 2 l«l h )J 

where 

A = ||| K (l + e- 2 l«l h ) 

1 + M 

Then with equations (A26) and the substitutions 

1/2 

cos|x=(^) J_i/ 2 (£ x ) 


(| x | > l) (A24b) 


(A25) 


(A26) 


3/2 

f(£) = - » A -(cosh 2£h + 1) 
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equations (A24) become 


J_ 1/2 (|x) d| (I x | < l) (A27a) 

and 

0 = f f(|)J_ 1/ / 2 (?x) d£ (|x[ > l) (A27b) 

Equation (A27b) is satisfied by the assumption 

OO 

f(?) = 3TiI { ' 1/2 I p m j 2m + l<5> 
m=0 


M 

3 + ii 




sinh 2£h + 2|h 


1 - M 

3 + Hu 


cosh 2£h + 1 


and further use of the method of reference 5 yields the familiar equations (A4), where 
now 


Lm,n - -2<2" * » t 


- 2 ?h (ir^) 


+ e 


-2|h 


cosh 2£h + 1 


J 2m+1 ( ^ J 2n+1^¥ 


The quantity of particular interest here is the moment My(x,0), given by 


M, 


' (x > 0) = -il 


" 2- n 

±S(i,0) - Ml 2 w(S,0) 
dy z 


cos £x di; 


With the appropriate substitutions, the moment is found to be 


ou 

My(x,0) = -M £ P m 




m 


m=0 L(x2-i) 1/2 [x + (x 2. l) 1/2 ] 


2m+l 


X 


00 1 - 2^h 




+ e 


-2£h 


cosh 2£h + 1 


J 2m+1^ cos d 9 (X>1) (A28) 


The infinite integral in equation (A28) can be proved convergent at x = 1, so that the 
ratio of the "moment intensity factors" is found to be 
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m 


which is equivalent to the ratio of the outer-fiber bending stress intensity factors, since 
the stress is given in terms of the moment by 



6My 

" 5 ~ 


Since the L m n values depend on Poisson’s ratio, so also does the stress inten- 
sity factor. It has been computed for four values of fi in the range 0 S ju = 0.5. The 
results are presented in figure 12 where, because of the relative insensitivity of Kjy/K^ 
to n, only the extreme values (for ix = 0 and g = 0.5) are shown. 
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PROBLEMS INVOLVING CONCENTRATED LOADS ON CRACK FACES 

This appendix contains most of the analytical details of the problems concerned with 
opening of a central crack in a strip by the application of equal and opposite concentrated 
loads. The problem involving free longitudinal edges has been given detailed treatment in 
the body of this report. For the sake of brevity, it will be referred to in this appendix, 
wherever possible, in the discussion of the remaining concentrated load problems. 

Crack Opening Mode - Concentrated Forces 

Zero normal displacement of longitudinal edges .- In this problem the crack faces 
are being separated by symmetrically applied concentrated forces of dimensionless mag- 
nitude P/Ea. The edges of the strip are free of shear stress and are restrained against 
normal displacement. The boundary conditions are 

v(x,h) = 0 
°xy( x >h) = 0 
°xy( x >®) = 0 


oy(x,0) = - 6(x) 

v(x,0) = 0 


In terms of <p, the transformed boundary conditions are 

= 0 
dy 3 

£%h) = 0 

dy 




(i*i<») 

(1*1 >l) 




> (Bl) 


oo 

6(x) = i j £ 2 0(£,O) cos £x d| 

(i*i< i) 

0 = j* cos £x 

J 0 dy 3 ^ 

(i*i>») 
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(B2) 


/ (B3) 


Substitution of equation (13) into the first three of equations (Bl) gives 

A = -g(e 2 l«l h + 2|«|h- l)'’ 

B = - me(e 2 l«th- x) l 

C = - £(2|||h+ 1 - e‘ 2 l®l h ) 

J 

with 

i = |l|(l-e- 2 l«l h ) 

Substitution of the identity (23) and equations (B2) into the last two of equations (Bl) 
yields 

!i(! ) l/2x ' i/2 «< x) = y; ora j.i/ 2<« o*!* 1 ) 

°.y; t ( 9J . 1/2 ( Wdl dx|>i) 

In this problem, 

f(!) = - ^(cosh 2£h - 1) 

G(t) - t/sinh 2|h + 2|h\ 

V cosh 2£h - 1 j 

In accordance with reference 5, it is assumed in equations (B3) that 

= fe * Z Pm J 2m+1^ 

m=0 


which satisfies the second of equations (B3) automatically. Further application of the 
method of reference 5 yields 

00 

P n + ^ Bjn^Pjn =1 (n = 0, 1, 2, . . .) (B4) 


m=0 


in which, as usual, 


L m ,n = 2(2n + 1) - l] J 2m+1 (l) J 2n+ i(S) ^ 


(B5) 
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After solution of equations (B4) for a sufficient number of the P n values, the stress of 
particular interest is found to be 


r 




1 + 2th - e-^ h 
cosh 2£h - 1 


J 2m+1^ cos d S/ 


(x > 1) 


from which the ratio of stress intensity factors is determined as 


K 

/=l + 2 X (-D m (Pm - 1) 


m=0 


(B6) 


The normal stress on a longitudinal edge is 


Oy(x,h) = 



sinh gh+ gh coshjh j (g) cos 
cosh 2£h - 1 2m+ l 


£xd£ 


It is symmetric in x. 

Clamped longitudinal edges .- In this problem, the edges of the strip are completely 
restrained against normal and tangential displacement. The boundary conditions are 

u(x,h) = 0 

v(x,h) = 0 

°xy(x,0) — 0 


<7y(X,0) = “ 6 ( X ) (l X I < !) 

v(x,0) = 0 (I x I > l) 

In terms of <j>, the transformed boundary conditions are 

2 - > 

^-4(4, h) + =0 

dy 2 

0(^,h) - (2 + n)b 2 ^(?,h) = 0 

* 0 1 

(Equations continued on next page) 
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IE = i J n £ 2 < MS>°) cos S x d£ (| x |< l) 

(M>1) 


0 = f ^-4(1,0) cos * 

J 0 dy 3 i 2 


j 


Substitution of equation (13) into the first three of equations (B7) gives 


A - - s(ft i e21 5 |h + 2 « 2| > 2 - 2 1« l h + TT7 + j 
B= - J ^(fr^ e2U|h - 2 i«i h+1 ) 

c = -5(r^ e ' 2U|h + 2 « 2h2 + 2 l { l h+ TT7r +1 ) 


(B7) 


(B8) 


where 


A = 



and 


K 


1 ~ M 

1 + jU 


By virtue of the identity (23) and equations (B8), the last two of equations (B7) take on the 
form of equations (B3), in which, for this problem, 


m = £ 


3 / 2 d 


sinh 2|h - 2|h 

1 + jU. 


GiQ = d ^ 


cosh 2|h + 2£ 2 h 2 + + l\ 

1 + a 1 + u 


sinh 2£h - 2£h 

1 + M 


/ 


Further use of the method of reference 5 yields equations (B4) and (B5). The stress of 
primary interest is 
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ou 

a y (x ’ 0) = ^ Z Pm 


Szl 1 


m 


■» W - i) 1/2 [ x ♦ u* - i) 1/2 l 


2m+l 


r 

J o 


00 1 + -y— r + 2|h + 2^ 2 h 2 + r— e" 2 ^ h 


1 + /i 


1 + M 


l-^Ji sinh 2£h - 2£h 
1 + /i 


J 2m+1^ ) cos d V (X > X) 


The ratio of stress intensity factors again is given by equation (B6). In this problem, 
the ratio of stress intensity factors depends on fj. because of the clamped- edge condi- 
tions. It has been computed for /i = 0 and p. = 0.5. 

The normal stress on an edge of the strip is 


o y (x,h) = - 


4P 


irEa(l + fi) 


2 



2 coshfti + (1 ± sinh jh J 2m+1 (£) cos |x d£ 

?-=-£ sinh 2£h - 2£h 
1 + fJi 


This stress is symmetric in x. 


52 



APPENDIX C 


CONVERGENCE PROOF FOR AN INFINITE INTEGRAL 


The integral to be investigated for convergence at x = 1 is 


lx=l 


(m = 0, 1, 2, . . .) 


1(1, h) = f 2 i 2 h 2 ± 2fo + 1 - e~ 2 ^ h . . 

Jo sinh 2£h + 2|h J 2m+1^ cos ?x d? 

where h is real and positive. Let the integral be written as 

r\ OO 

= J o H (^ h ) J 2m+1^ cos * (Cl) 

The function H(£,h) is a positive steadily decreasing function, for all h > 0, in the 
interval 0 < £ < °o. it is well-known that 


and 


so that 


2m+l 


cos £ 


(£) 


< 1 


i 1 


(m = 0, 1, 2, . . .) (C2) 


(C3) 


J 2m+1® cos «| 


< 1 


(C4) 


Assume that the integral 


J H(^,h) d£ 


exists. If the integral (C5) exists, so does the integral 


r\ OO 

J 0 H(«,h) |j 2m+1 ({) co S { 


d£ 


(C5) 


(C6) 


which follows from a simple comparison test for infinite integrals (see ref. 9, p. 71), by 
virtue of the inequality (C4). By use of a theorem on absolute convergence of infinite 
integrals (see ref. 9), it further follows that the integral 


r "(«’») ■wi«> c ° s « 


53 



APPENDIX C 


exists, if the integral (C6) exists. Therefore, all that remains for proof of convergence 
of equation (Cl) is to prove the convergence of integral (C5). 

Now H(£,h) can be written as 


Note that H(£,h) = 1 at £ = 0. For £ > l x , where ^ is the finite (real) root of 

(4{h-e‘ 2 « h )e- 2 « h = 0 


it is seen that 

H(£,h) < 2e" 2 ^ h (2^ 2 h 2 + 2£h + 1 - e" 2 ^ 1 ) 

It is then permissible to write, for e > 0, 

f H(£,h) d£ < f 1 H($,h) d£ + 2 j* e 2 ^ h (2| 2 h 2 + 2£h + 1 - e" 2 ^ h ) d£ ^ 


Both integrals in the right-hand side of the inequality (C7) are convergent. The first 
integral is convergent because H(|,h) contains no singularities, integrable or other- 
wise, in the range of integration, and the second integral converges because every term 
in the integrand is of negative exponential order. Therefore, the integral (C5) exists, 
and the proof is complete. 
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